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ABSTRACT 


The objective of this research was to provide the Office of Naval 
Research (ONR) with an independent evaluation of the structural adequacy of 
Lockheed's SLICE hull design. SLICE is an offshoot of the SWATH ship 
technology, with four lower hulls instead of the SWATH's twin hulls. An advanced 
finite element code was used in the dynamic and static analysis of the SLICE's 
unique hull form. This thesis demonstrates the practical value of utilizing finite 
element modeling as a tool in modern ship design. 
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I. INTRODUCTION 


A. BACKGROUND 


Design of naval vehicles has long been an art passed down from generation to 
generation. Structural design rules of the various classification societies have evolved 
through years of operational experience. Because of the complexity of hydrodynamic 
loads, hull structural designs have relied heavily upon static analysis as modified by 
empirical data to account for actual dynamic conditions. Although this traditional 
approach to hull design has proved to be adequate, it does not efficiently utilize material 
consistent with the loads and failure modes hulls need to resist. With the quantum 
increase in computing power made available during recent years, there is no longer any 
reason not to overhaul these traditional methods and to incorporate the latest technology 
into new ship designs. 

The basic hull structure of a ship must be designed to resist yielding, fatigue and 
buckling throughout its service life. A finite element analysis of the structure using 
expected hydrodynamic loading conditions allows the designer the flexibility to 
effectively vary structural member sizes (scantlings) throughout the ship, based upon the 
resultant stresses and stress concentrations. This dynamic approach allows an optimum 
structural weight distribution and minimizes the opportunity for structural failure. The 
limitation to such an approach is the current lack of available data to accurately predict 
sea wave spectra throughout the life of a ship; however, this same limitation also plagues 
the traditional design approach. 

A lack of historical data makes the finite element model approach especially 
appropriate when considering radically different hull forms as opposed to the 
conventional monohull. This requirement has become increasingly evident over the past 
several decades with the growing use of the Small Waterplane Area Twin Hull (SWATH) 
ship design. This innovative hull form requires new hydrodynamic pressure load criteria 


since the existing methods developed for monohull ships do not apply. Unlike the 


limiting longitudinal loading condition which causes hogging and sagging in monohuils, 
the primary wave loads acting on the SWATH’s submerged hulls and water piercing 
struts are in the transverse direction. The remainder of this paper demonstrates these 
principles with a novel ship design known as the SLICE, a variant of the SWATH 


design. 


B. DESCRIPTION OF THE SLICE ADVANCED TECHNOLOGY 
DEMONSTRATOR (ATD) 


The SLICE Advanced Technology Demonstrator (ATD) is a unique surface ship 
hull form designed by Lockheed Missiles and Space Company (Marine Systems Division) 
and sponsored by the Office of Naval Research (ONR). The design concept for SLICE is 
an offshoot of the Small Waterplane Area Twin Hull (SWATH) ship technology 
developed during the past three decades. Its unique hull form consists of four submerged 
“torpedo-like” lower hulls or pods connected by means of struts to the main cross 
structure located above the water’s surface. The forward pair of eight foot diameter hulls 
is inset by three feet with respect to the aft pair. Figure | depicts the outboard profiles 
and principal dimensions of the SLICE ATD. 

The primary advantage of the SLICE/SWATH over conventional monohulls is its 
enhanced seakeeping ability, this is especially true in heavy seas. Buoyancy is 
predominantly provided by the fore and aft set of hulls which remain well below the 
water’s surface, effectively isolating the ship from the dynamic buoyancy effects of 
waves. The ship's four water piercing struts present a much smaller waterplane area than 
a conventional monohull of similar displacement which make it less sensitive to wave 
actions. Other advantages include [Ref. 2]: 

e Reduced deck wetness 

© Reduced wave slamming 

© Increased crew effectiveness and safety in high sea states 


e Ability to maintain speed in high sea states 
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Figure 1. Outboard Profile and Principal Dimensions of the SLICE ATD, from Ref. [1] 


Despite its superior seakeeping ability, the SWATH ship is not without 
limitations. Because of its small waterplane area, the SWATH displays hydrostatic 
properties markedly different from a monohull’s. Its typical tons-per-inch immersion 
(TPD) is 15-30 percent that of a similarly displaced monohull’s TPI. Similarly, its 
moment-to-change trim one inch (MTI) is 10-20 percent that of a monohull’s. 
Consequently, due to its increased sensitivity to weights and moments, the SWATH’s 
displacement and weight distribution becomes an important factor during the preliminary 


design process. [Ref. 3] 


C. HYDRODYNAMIC LOADS ON SWATH SHIPS 
1. Primary Loads 


The structural configuration of SWATH ships is strongly influenced by the 
hydrodynamic forces of waves. The passage of waves creates a variable pressure 
distribution that is generally different on its twin hulls (or the four hulls of a SLICE). A 
typical pressure distribution resulting from beam seas is shown in Figure 2. The resultant 
forces from these waves are maximum in beam seas, approximately 50-60 percent of 
maximum in bow or quartering seas, and small in head and following seas. These forces 
which alternately pry the hulls apart and squeeze them together are the primary hull girder 
loads for SWATH ships. The magnitude of the alternating transverse side forces 
generated varies with the seas encountered and hull form selected; however, 
measurements and model tests have shown that they are comparable in size to the 
displacement of the ship. [Ref. 3] 

Due to a phenomenon known as shear lag, the transverse bulkheads and adjacent 
shell plating must absorb most of the shear and bending moments caused by the primary 
loads. Since the shell plating between the bulkheads is much more flexible than the 


transverse bulkheads in the athwartships direction, it provides little resistance to the 


wave induced forces. Consequently, the placement and configuration of transverse 


bulkheads becomes a critical factor in the design of SWATH ships. (Ref. 3] 




















Figure 2. Primary Load Pressure Distribution in Beam Seas, from Ref. [3] 


2. Secondary Loads 


The principal types of secondary loads a SWATH ship encounters are wet deck 
slamming and strut side slamming. Slam pressure loads are local in nature and occur at 
all locations on the wet deck; however, they are generally most severe at the ends due to 
the ship’s pitch motion. Prediction of strut side slamming is difficult due to complexity 
of flow near the hull. Very little data has been collected to date concerning slam 


pressures on SWATH ships which make this a difficult area to analyze. [Ref. 3] 


D. SCOPE OF RESEARCH 


This thesis documents the method used in the development of the finite element 
based representation of the SLICE ATD hull. An advanced finite element code called 
LDEAS (Integrated Design Engineering Analysis Software) developed by Structural 
Dynamics Research Corporation is utilized as a design tool for the analysis. The software 
is used to solve for the normal modes of oscillation as well as the linear static response of 
the SLICE ATD finite element model. Subsequent chapters will discuss the theory 
behind finite element modeling, model development and design decisions, analysis 


results, conclusions and suggestions for further research. 


IL. FINITE ELEMENT MODELING 
A. BASIC DESCRIPTION 


Finite Element Analysis (FEA) is a numerical process used to approximate 
displacements and stresses within a structure. The method involves dividing a complex 
structure into a grid of discrete elements. Each element consists of a simple polygon for 
which the finite element program assembles a set of governing equations in the form of 
stiffness and mass matrices. The unknowns for each element consist of the displacements 
at the “nodes” or points at which the elements are connected together. There will be one 
equilibrium equation written for each nodal degree of freedom; Equation (1) illustrates 
the equations of motion for a dynamic system in matrix form. Given the known forces, 
structure stiffness, structure mass, viscous damping and boundary conditions, the 
equations of motion can be solved to yield the displacement at the nodes. The stresses in 


each element can then be calculated from the nodal displacements. [Ref 4] 


(M1(q} +[C](4} +(K]{q} = (F} (dy 


where [M]= mass matrix, [C]= viscous damping matrix, [K]= stiffness matrix, 
{q}= generalized nodal displacement vector, {q}= generalized nodal velocity vector, 
(G} = generalized nodal acceleration vector, and {F}= nodal force vector. 

Finite Element Modeling (FEM) is the process of creating a structure through the 
use of these nodes and elements. The finite element model is a mathematical 
representation of a structure and may not necessarily look like the actual structure. 
Different types of elements are used to model the various parts of the structure including 
thin shell, beam, solid, rigid and lumped mass elements. The accuracy of the results 
will depend upon the proper selection and combination of these element types. 
Additionally, the assumptions made for loads, boundary conditions, and material 


properties significantly affect the accuracy of the solution. In general, the solution will 


be more accurate as the model is subdivided into smaller elements, however, solution 
accuracy does not come without a price. The price you must pay comes in the form of 
increased computer processing time and data storage requirements which often don’t 
justify the need for an incremental increase in accuracy. The art of finite element 
modeling is, therefore, to discretize the structure into an appropriate number of finite 
elements to provide adequate accuracy without using an excessive number of elements. It 
is up to the modeler to apply sound engineering judgement when selecting element types, 
material properties, loads, and boundary conditions in order to gain useful insight into 


the response of a complex structure. 


B. OVERVIEW OF FINITE ELEMENT MODELING USING I-DEAS™ 


I-DEAS™ (Integrated Design Engineering Analysis Software) is an integrated 
package of mechanical engineering software tools used for solid modeling, drafting, and 
finite element simulation. FEM using I-DEAS is an iterative process with three main 
stages [Ref 5]: 


©  Pre-processing includes developing the model geometry, assigning physical 
and material properties, applying loads and boundary conditions, and 
checking the model quality. 


« Solving the model is done in the I-DEAS Model Solution Task. I-DEAS 
Model Solution can perform linear and non-linear statics, dynamics, and 
buckling analysis. The model information can also be exported to several 
popular external finite element solvers for analysis. 


© Post-processing involves plotting deflections and stresses, and comparing the 
results to predetermined failure criteria. I-DEAS shows the results in different 
display formats- deformed geometry, contour plots, and others- to give 
insight into how to improve the design. Post-processing also allows the 
engineer to check for errors which may not have been detected during the 
model building process. 


Figure (3) represents the sequential steps required to complete the FEM process using 
I-DEAS. Details on the pre-processing stage are included in Chapter III. Details on 


solving and post-processing the model are discussed in the following section. 
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Figure 3. The FEM Process, after Ref [5] 





C. MODEL SOLUTION 


The I-DEAS Model Solution Task is an integrated finite element structural 
analysis solver which can be used for both dynamic and linear static analysis. Both of 
these analysis techniques provide valuable insight into the response of complex 


structures. 


1. Normal Mode Analysis 


Normal Mode analysis is a technique for predicting the undamped natural 
frequencies and mode shapes of vibration for a structure. Mode shapes and natural 
frequencies can be used to identify transient load points and frequencies which may 
generate undesirably large structural responses to dynamic inputs. The response of 
structures can often be assumed to be a combination of the mode shapes corresponding to 
the lowest natural frequencies. [Ref 6] 

In the absence of damping and external forces for a multi-degree-of-freedom 


system, the equations of motion reduce to: 


(M1{4} + [K1{q} = (0} (2) 


The solution to this set of second-order differential equations is not an easy task, even 
when the equations are linear with constant coefficients. This is because the equations 
are dependent upon one another through coupling terms and must be solved 
simultaneously. The equations of motion are typically coupled through stiffness or mass 
terms, known as elastic and inertial coupling, respectively. By transforming the physical 
coordinate system to one which diagonalizes the stiffness and mass matrices 
simultaneously, the equations of motion can be uncoupled both elastically and inertially. 
Such coordinate systems do exist and are known as natural coordinates. The process of 


diagonolizing the matrices is sometimes referred to as modal decomposition. Once the 


diagonalized matrices are obtained, the solution becomes one of solving for n 
independent homogeneous equations, where n is the number of degrees of 
freedom (DOF) for the system. The general motion of an undamped linear system can 
then be represented as a linear combination of its natural modes of vibration. [Ref 7] 

The I-DEAS Model Solution Task has three different algorithms available to 
solve for natural frequencies and mode shapes. These are: Simultaneous Vector 
Iteration (SVI), Guyan Reduction, and the Lanczos method. All three methods use some 
form of modal decomposition to obtain the results. The SVI method was chosen for use 


in this thesis due to its efficiency when working with a large DOF system. (Ref 5] 


2. Linear Static Analysis 


The more traditional approach to finite element analysis is through the use of 
linear statics. Linear statics is a technique used to predict the displacement of a structure 
when a static load is applied within the elastic range of the structure's material properties. 
It uses the principle of virtual work to equate the virtual work of external forces and 
moments to the virtual work of internal stresses. Once the nodal displacement of the 
finite element model is known, other important parameters such as stress can be 


calculated. The equations of motion in matrix form for linear statics are simply: 


[A1{9 = (A (3) 


Again, the number of equations corresponds to the number of DOF for the 
system. [Ref 6] 

To obtain reasonable results for a linear static analysis, proper boundary 
conditions must be applied to the structure. A boundary condition set consists of both 
restraints and loads. Improperly restraining a structure can lead to misleading and 
erroneous results. A model must be restrained against all six rigid body modes (three 


translations and three rotations). If rigid body motion is a possibility in a structure (as in 


a ship or other vehicle), it ast be restrained using kinematic degrees of freedom. 
Kinematic degrees of freedom act to restrain a vehicle in space when no physical 
restraints exist. [Ref 4] 

Finally, the results of the normal mode analysis can be used to help identify 
probable limiting load cases for use in a linear static analysis. External forces may be 
obtained from existing data bases such as data obtained from model or full scale tests of 


similar structures. 


D. POST-PROCESSING 


Post-processing is the visual display and interpretation of data from the model 
solution stage. Two of the most valuable displays available within I-DEAS are the 
contour plot and the deformed geometry plot. Contour plots display the distribution of 
results using contour lines and bands of color representing the data. This is especially 
useful when evaluating the stress distribution in a structure. Deformed geometry plots 
display the model in its deformed condition, i.e., as if the displacements were bending 
the model. Because most deformations are small compared to the size of the model, the 
displayed deformations are magnified by a scaling factor to allow for meaningful 
interpretation of data. [Ref 5] 

Another useful feature of I-DEAS is its ability to animate the results. Animation 
gives the illusion of motion to the deformed geometry by succesively stepping through a 
series of frames, starting with the undeformed geometry and ending with the 100% 
deformed configuration. This provides valuable insight when evaluating the mode shapes 


of a structure. [Ref 5] 


Il. MODEL DEVELOPMENT 


The purpose of this chapter is to describe the development of the SLICE ATD 
finite element model used for this study. The important design decisions and engineering 
simplifications used during construction of the model will be discussed. The finite 
element model of the SLICE hull was created by using Lockheed’s detailed design 
drawings as a reference. Figure 4 shows the profile and main deck view from these 
drawings. Figure 5 shows similar views of the completed SLICE finite element model for 


comparison. An isometric view of the completed model is illustrated in Figure 6. 


A. PHYSICAL AND MATERIAL PROPERTIES 


Before creating elements, the physical and material properties of the elements 
must be defined. Physical properties are the geometric properties of an element. The 
only physical property varied for this model was the element thickness, which varied 
from 0.09 inches for the non-tight (NT) floors to 0.75 inches for the main transverse 
beams. 

Material properties are those properties which are unique to a specific material. 
They include Young’s Modulus, yield strength, Poisson's ratio, and material density. 
The hull structure of the SLICE is to be constructed from 5083-H32 Aluminum for plates 
and 5083-H112 Aluminum for shapes. Material properties are listed below in Table 1. 




















Material Type 5083-H32 5083-H112 
Elastic Modulus 10.3 x 10° psi 10.3 x 10° psi 
Yield Strength (min) 31,000 psi 24,000 psi 
Poisson’s Ratio 0.33 0.33 
Density 096 Ibmivin® 096 lbm/in® 














Table 1. Material Properties of 5083 Aluminum, from Ref [8] 
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Figure 4. Structural Arrangement of the SLICE, from Ref [9] 
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Figure 5. Finite Element Model of the SLICE 
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Figure 6. Isometric View of the SLICE Model 





B. MESH GENERATION OF THE SLICE 
1. Node and Element Creation 


Each node is a point located in 3D space where elements are connected, loads are 
applied, and restraints are imposed. Each node has up to six degrees of freedom, 
depending on the type of element used. Nodes are located such that the geometry of the 
model approximates that of the actual structure. Nodes for the SLICE model were 
generated manually by keying in their coordinates or by copying and reflecting existing 
nodes. Each transverse frame (see Figure 4) of the SLICE was modeled as a set of 
nodes. [Ref 4] 

Prior to generating any elements, the decision was made to construct the basic 
hull form from thin shell linear quadrilateral and triangular elements. Thin shell elements 
can be used to accurately represent a structure whose thickness is small with respect to its 
other dimensions. No one element type is best for all situations; however, thin shell 
elements most closely model the response of stiffened hull plating. Although cubic and 
parabolic elements provide more accuracy than linear elements, the small gain in 
solution accuracy was not deemed necessary for this study. Linear elements have two 
nodes along each edge, parabolic have three, and cubic have four. [Ref 4] 

Once the nodes are generated, the elements can be created. Each element is 
created by sequentially selecting a set of three or four nodes; i.e., three nodes for a 
triangle or four nodes for a quadrilateral. Adjacent elements share a pair of common 
nodes; thereby forming a “grid” or “mesh” of elements representing the components of 
the actual structure. 

Once the basic hull form was completed, lumped masses were added to the model 
to account for nonstructural mass components not represented by the thin shell elements. 
These include the propulsion plant, the payload structure, the deck house, the auxiliary 


machinery, and fuel oil. The Appendix contains a summary of the SLICE ATD weight 


distribution. The lumped masses were connected to the thin shell elements through a 





network of rigid elements. Rigid elements are massless, infinitely stiff members used to 
distribute the mass loading of the lumped masses. 

In summary, the SLICE ATD finite element model contains 10,262 nodes, 
11,562 linear thin shell elements, 456 rigid elements, and 23 lumped mass elements. 


The model has 58,875 degrees of freedom. 
2. Model Components 


a. Transverse Bulkheads 


The transverse bulkheads provide the main structural strength for the 
SLICE. The SLICE contains multiple transverse bulkheads; however, only two 
bulkheads provide its primary load carrying capacity. These primary bulkheads are 
referred to as the forward and aft spars which are located at frames 6 and 35, 
respectively (see Figure 4). The SLICE incorporates a unique design in the construction 
of these two bulkheads, henceforth referred to as the forward and aft spars. Each spar 
employs use of a tapered I-beam which is curved by 90 degrees as shown in Figure 5. 
The flange thickness of the curved beam is 0.75 inches. Its web thickness varies from 
0.50 inches at the curve or haunch area to 0.25 inches at the ship centerline and strut. 
The longitudinal girder shown in Figure 4 has a flange and web thickness of 0.25 inches. 


The model’s thin shell element thickness was varied to reflect these dimensions. 


6. Stiffened Hull Plating 


The hull plating was modeled using 0.25 inch thin shell elements. The plating is 
strengthened through the use of both transverse and longitudinal stiffeners. To minimize 
the degrees of freedom of the model, the smaller stiffeners’ cross sections (i.e., those 


with a depth of 6 inches or less) were idealized as a wedge shape vice the typical Tee 
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Figure 7, SLICE’s Primary Transverse Bulkheads 








shape. This was done by matching the moment of inertia (in the direction of bending) 
and the center of gravity of the wedge to that of the corresponding Tee. Figure (8) 


illustrates this concept. A typical section of stiffened hull plating is shown in Figure (9). 
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Figure 8. Idealization of Tee Shaped Stiffeners 
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Figure 9. Stiffened Hull Plating of Forward Main Deck 


C. MESH QUALITY CHECKS 


I-DEAS contains several mesh quality checks to help identify modeling errors. 
Typical errors which can be checked include element connectivity problems, highly 
distorted elements, and duplicate nodes or elements. [Ref 4] 

The free edge check plots the free edge of elements which are not connected to 
other elements. Normally, this plots only the outer boundary of the model; however, if. 
elements are not properly joined, an extra line will appear when conducting this check. 
This represents a “crack” or flaw in the model. [Ref 4] 

A second useful check is the element distortion check. In finite element 


modeling, a “perfect” element is represented by a square or equilateral triangle. Any 





deviation from this ideal condition results in a distorted element. Highly distorted 
elements should be avoided, especially in important areas such as high stress 
locations. [Ref 4] 

Finally, I-DEAS can check for coincident nodes and elements. This feature is 
especially important when building large or complex models. Having duplicate nodes or 


elements can produce erroneous data or cause the finite element program to fail. 


D. BOUNDARY CONDITIONS 


Prior to solving the finite element model, boundary conditions must be applied. 
Boundary conditions can consist of both restraints and structural loads. For the Normal 
Mode Analysis, only restraints are used. Even if a model has rigid body motion, it must 
be held in space by restraints so that it is not free to move in any direction. Restraints are 
applied at individual nodes and can have up to six values at each node, three translations 
and three rotations. In the case of the SLICE, six rigid body modes are possible; 
therefore, the model must be immobilized by a minimum of six restraints referred to as 
kinematic degrees of freedom. A total of nine translations were selected at three separate 
nodes to ensure the model was sufficiently restrained. (Ref 4] 

For the Linear Static Analysis, structural loads are required in addition to the 
restraints. For the SLICE model, structural loads include gravitational forces, buoyancy 
forces and hydrodynamic loads due to waves. Gravitational forces are distributed 
throughout the model and vary depending upon the individual mass of each element. 
Buoyancy forces are distributed around each of the four submerged hulls at nodes which 
are close to the ship’s center of buoyancy. 

By far the most complex forces to simulate are those due to the sea. One useful 
method available to predict the maximum lifetime load for SWATH ships is presented in 
references [11] and [12]. This method develops a design side load algorithm which 


predicts the single maximum side force a SWATH ship would experience during a 


2 


20 year lifetime in the North Atlantic Ocean. A summary of the algorithm developed in 
references [11] and [12] is presented here. The units of displacement are Ltons and the 
lengths are in normalized units of feet divided by the cube root of 


displacement (Feet/Ltons'”) 
F=DTLA (4) 


where F = Maximum Lifetime Side Load (Ltons) 

A= Displacement (Ltons) 

D= Size Factor = 4.474 - 0.913*Log(A) 

L = Length Factor = 0.75 + 0.35*tanh(0.5*Ls - 6) 

where Ls = Strut Length (Feet/Lton') 
[for tandem strut vessels, Ls is the sum of the forward and aft strut lengths] 

T = Draft Factor = 0.5319*Draft (Fee/Lton'*) 

Inputting the appropriate SLICE parameters into Equation (4) yields a maximum 
lifetime side force of approximately 280 Ltons. This side load is overestimated for small 
SWATH ships that are restricted to coastal operations because of the assumption of a 
maximum wave height of 50 feet used to derive the algorithm [Ref 3]. A reduced initial 
estimate of a 200 Lton side force was used for this study; however, further analysis 
needs to be performed to refine this value. 

Another finite element study determined that the variations in side load 
distribution had relatively little effect upon the stresses in SWATH structures [Ref 13]. 
Variations in the loading redistributed themselves throughout the structure [Ref 13]. The 
distribution of side forces and buoyancy forces chosen for the current study is illustrated 
using the aft starboard hull in Figure 10. A trapezoidal side force distribution was 


assumed with the resultant force centered about the mid-draft level. 
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Figure 10. Typical Distribution of Side Forces and Buoyancy Forces 
(Note: Arrows represent direction and not magnitude) 





IV. DISCUSSION OF RESULTS 
A. NORMAL MODE ANALYSIS 


The SLICE model's first four flexible modes were solved using the Simultaneous 
Vector Iteration (SVI) method of I-DEAS’ Model Solution Task. The mode shapes were 
evaluated in the Post-Processing Task using I-DEAS’ animation feature. A_ brief 


description of each is provided here. 


1, First Mode 


The first mode occurs at a natural frequency of 4.7 Hz and exhibits the classical 
squeezing and prying motion expected of SWATH ships. The forward and aft pods are in 
phase with each other throughout the cycle. The forward pods also exhibit a torsional 
response during the cycle. Figure 11 shows the deformed geometry plot of this mode 


shape during the pry cycle. 


2. Second Mode 


The second mode occurs at a natural frequency of 5.1 Hz with the aft pods 
displaying a squeezing and prying motion. The forward pods exhibit a “racking” motion 
in which the two pods travel in the same direction. The forward pods travel to starboard 
during the aft pod’s pry cycle and travel to port during the squeeze cycle. Once again, a 
twisting motion is observed on the forward pods. Figure 12 shows this mode shape 


during the aft pod’s pry cycle. 
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3. Third Mode 


The third mode is very similar to the second mode and occurs at a natural 
frequency of 5.5 Hz. The main difference between the two modes is that the racking 
motion of the forward pods moves to port during the pry cycle and to starboard during the 


squeeze cycle. Figure 13 displays this mode shape during the pry cycle 


4. Fourth Mode 


The fourth mode has a natural frequency of 6.1 Hz and displays the classical 
longitudinal bending response of a monohull. There is no pod torsion as observed with 


the previous modes. Figure 14 illustrates this mode shape. 


B. LINEAR STATIC ANALYSIS 


A linear static analysis was conducted on the model based upon the results 
obtained during the normal mode analysis. Three load cases were considered using a 
200 Lton side load superimposed with a 168 Lton buoyancy force and a 168 Lton 
gravitational force. The 200 Lton side force was divided equally between the forward and 
aft pods. The three load cases evaluated were: 

(1) Prying force applied to the forward and aft pods, 

(2) Squeezing force applied to forward and aft pods, and 

(3) Racking force applied to the forward and prying force applied to the aft pods. 


1. Load Case One (Prying Force) 


Figure 15 displays the Von Mises stress contour plot for the SLICE’s external 


skin. The units of stress are expressed in pounds force per square foot (Ibf/ft?) on the 
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Figure 11. Deformed Geometry Plot of First Mode (4.7 Hz) 





































































































Figure 12. Deformed Geometry Plot of Second Mode (5.1 Hz) 
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Figure 13. Deformed Geometry Plot of Third Mode (5.5 Hz) 
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Figure 14. Deformed Geometry Plot of Fourth Mode (6.1 Hz) 
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color bar. Due to shear lag, the transverse bulkheads and adjacent shell plating absorb 
most of the shear and bending moments caused by the loads. The stress contour plot for 
the tranverse bulkheads is shown in Figure 16. The forward and aft spars were designed 
to carry most of the primary load; and that is in fact what occurs. The forward and aft 
spar stress contour plots are presented in Figures 17 and 18, respectively. It is interesting 
to note that the high stress areas are predominately in the struts and not in the haunch 
area (as in previous SWATH designs). This can be attributed to the unique design of the 
spars. The curved I-beam reduces stress concentrations and is also strengthened by using 
a 0.50 inch.web thickness in the haunch area. Because the struts only have a 0.25 inch 
web thickness, the stress is higher in that region. The stress in the strut region (=21 ksi) 


comes dangerously close to the matrerial’s yield stress (24 ksi). 


2. Load Case Two (Squeezing Force) 


The results for the squeeze cycle are very similar to the pry cycle. The Von Mises 
stress contour plots are presented in Figures 19 through 22. Once again, the struts of the 
forward and aft spars are the limiting stress regions with a maximum stress of 


approximately 21 ksi. 


3. Load Case Three (Racking Force Forward/Prying Force Aft) 


This unique loading condition places a torsional moment along the SLICE’s 
starboard longitudinal girder. As shown in Figures 23 and 24, the aft starboard pod and 
adjacent transverse bulkheads absorb most of the shear and bending moments. As 
observed from Figures 25 and 26, the stress in the aft spar is approximately twice the 
magnitude of the forward spar. The limiting stress for this loading condition is in the aft 


spar’s strut with a maximum value of approximately 24 ksi. 
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Figure 15. Von Mises Stress Contour Plot of Entire Ship (Prying Force) 
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Figure 16. Von Mises Stress Contour Plot of Transverse Bulkheads (Prying Force) 

















Figure 17. Von Mises Stress Contour Plot of Forward Spar (Prying Force) 











Figure 18, Von Mises Stress Contour Plot of Aft Spar (Prying Force) 
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Figure 19. Von Mises Stress Contour Plot of Entire Ship (Squeezing Force) 
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Figure 20. Von Mises Stress Contour Plot of Transverse Bulkheads (Squeezing Force) 











Figure 21. Von Mises Stress Contour Plot of Forward Spar (Squeezing Force) 
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Figure 22, Von Mises Stress Contour Plot of Aft Spar (Squeezing Force) 
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Figure 23. Von Mises Stress Contour Plot of Entire Ship 
(Racking Force Forward/Prying Force) 
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Figure 24. Von Mises Stress Contour Plot of Transverse Bulkheads 
(Racking Force Forward/Prying Force Aft) 














Figure 25. Von Mises Stress Contour Plot of Forward Spar 
(Racking Force Forward/Prying Force Aft) 









Figure 26. Von Mises Stress Contour Plot of Aft Spar 
(Racking Force Forward/Prying Force Aft) 
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V. RECOMMENDATIONS AND CONCLUSIONS 


A. RECOMMENDATIONS 


Based upon the results of the normal mode and linear static analysis of the SLICE 
finite element model, the following recommendations and suggestions for continued 
research are presented: 

« The SLICE prototype’s lower strut web thickness should be increased to 

0.50 inches. This will provide an increased margin of safety with minimal 
weight increase. 

e The analysis showed that the primary loads were effectively transmitted to the 
forward and aft spar as intended. The finite element model should be refined 
to provide a finer mesh in these high stress areas prior to conducting a more 
detailed analysis. 

e Further research is required to define realistic hydrodynamic loading criteria. 
This is required for both primary and secondary loads. 

e Using the expected hydrodymamic forces, evaluate the dynamic stresses 
imposed upon the SLICE structure with a solver such as the I-DEAS Forced 
Model Response function. 

e Finally, when operational, instrument the SLICE prototype vessel to validate 


the results of the finite element model. 


B. CONCLUSIONS 


Structural design of SWATH/SLICE ships is a complex, iterative process. The 
intent of this thesis was to demonstrate the value of using a finite element model as a tool 
in the design process. In the process of conducting the thesis, it was discovered that the 


most difficult and often confusing part of finite element modeling is the proper 
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application of boundary conditions. No matter how well one models a structure, without 
accurate representation of forces and restraints, the results are not reliable. With the 
processing capability of today’s computers, there is no reason to oversimplify these 
sometimes complex forces. Therefore, ship designers should seek to conduct further 
research into applying various realistic dynamic loading conditions which a ship may 
experience throughout its lifetime. Only then will designers be able to thoroughly and 


accurately evaluate design alternatives. 
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WEIGHT STATEMENT 


Summary 
Departure Cond 
ATO 
















































































Desenption TOTAL WEIGHT 
(los) «ny (ft) (ft) (tt) 
100 |Hull Structure 125183 33.26 | 0.03 | 18.76 
200 | Propulsion Plant 48228 22.73 | -0.05 | 9.95 
300 | Electric Plant 12651 65.99 | 5.18 | 21.52 
400 Instrumentation & Control 2144 32.37 | 0.00 | 30.46 
500 [Auxiliary Systems 24981 42.78 | -298 | 15.26 
600 _|Outfit and Furnishings 12439 25.12 | 0.22 | 24.58 
LIGHTSHIP 225627 33.44 0.02 17.07 
M10 |Margin - Weight (10.0%) 22563 10.07 | 33.44 | -0.02 | 17.07 
Margin - KG (6.7%) 1.14 
LIGHTSHIP w/MARGIN 248189 | 110.80 | 33.44 | -0.02 |) 18.22 
Fit |Crew (5) 800 0.36 20.00 | 0.00 | 41.75 
F41 | Fuel Oil (440 NM @ 30 Ktsi) 51758 23.11 | 43.38 | 0.00 | 10.53 
F46 — |Lube Oil (10 gal) TE 0.03 | 38.00 | 0.00 | 21.00 
F51 | TrinvBallast Water - Total o 0.00 8.20 0.00 3.73 
F52 | Potable Water (400 gal) 3328 1.49 | 25.00 | 0.00 | 27.00 
F52 [Payload gag98 42.41 | 51.07 | 0.04 | 92.19 
FULL LOAD CONDITION 399151 
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